Exact Mappings between Classical and Topological Orders in Two Dimensional Spin 

Models 
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Motivated by the duality between site-centered spin and bond-centered spin in one-dimensional 
system, which connects two different constructions of fermions from the same set of Majorana 
fermions, we show that two-dimensional models with topological orders can be constructed from 
certain well-known models with classical orders characterized by symmetry-breaking. Topology- 
dependent ground state degeneracy, vanishing two-point correlation functions, and unpaired Majo- 
rana fermions on boundaries emerge naturally from such construction. The approach opens a new 
way to construct and characterize topological orders. 
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Recently, topological orders have attracted intensive 
interests for different reasonsi^i^iiii^i^ii. The best studied 
example of topological order is the fractional quantum 
hall (FQH) states^. All different FQH states have the 
same symmetry. Unlike classically ordered state, FQH 
liquids cannot be described by Landau's theory of sym- 
metry breaking and the related order parameters^jl. A 
new theory of topological order is proposed to describe 
FQH liquids^. New nonlocal quantities, instead of local 
order parameters, such as ground state degeneracj^, the 
non-Abelian Berry's phasc"^ and topological entropy^ilS^ 
were introduced to characterize different topological or- 
ders. Topological ordered systems have also been de- 
signed and studied in the context of quantum computa- 
tion as an realization of potentially fault-tolerant quan- 
tum memory and quantum computationSi^iii. It is the 
non-locality of the topological orders that significantly 
reduces the effect of de-coherence^^. 

Theoretically, a number of soluble or quasi-soluble 
models which capture the topological orders have been 
proposed and studied^'^'^ d^d'^dSd^ . However, unlike the 
conventional orders which are entirely characterized by 
broken symmetries, the topological orders have not been 
characterized in a universal way. In fact, topological or- 
ders have to be studied case by case in different mod- 
els. Recently, it has also been pointed out that the spec- 
trum is completely inconsequential to topological quan- 
tum order^^ and hidden order parameter has been sug- 
gested in Kitaev model on honeycomb latticei^. In this 
work, we show a new way to characterize topological or- 
ders, which is based on well-known conventional models. 
First, we show that it is possible to map a model with 
topological order to a model with a local order param- 
eter in certain physical realizations through a non-local 
duality transformation. A local order parameter descrip- 
tion of topologically ordered systems is potentially use- 
ful. For instance, thermodynamic properties and energy 
spectrum can be easily computed in terms of classical or- 
der parameters. Second, we show that topologically or- 
dered systems can be constructed or designed from well- 
studied classically ordered states by including a topolog- 



ical boundary term associated with the lattice topology. 
In such a construction, topological properties are mani- 
festly presented. The result provides a novel approach for 
easier and/or better design of physical implementations 
of topological orders for quantum computation, starting 
from ordered systems well-understood in the framework 
of Landau's symmetry breaking theory. 

We start with examining a well-known nonlocal trans- 
formation, namely, the duality between site-centered 
spin and bond-centered spin in one-dimensional spin-1/2 



systemic, 
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The spin operators ti on the original lattice can be 
fermionized by a Jordan- Wigner transformation 
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where A{n) and B{n) are Majorana fermions on site n. 
Fermions can be defined as c(n) = [A{n) + iB{n)]/2. The 
duality transformation of Eq.(IT]) and Eq.([2]) now reads 



Hzin) ^ iB{n)A{n + 1) 
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B{n), (7) 



which is another Jordan- Wigner transformation if we in- 
troduce a new set of fermions d{n) — [B{n) + iA{n + l)]/2 
on the dual lattice. It is thus clear that the duality 
transformation connects two different constructions of 
fermions from the same set of Majorana fermions, as il- 
lustrated in FigU The duality now appears as a very lo- 
cal transformation. However, in terms of spin operators. 
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FIG. 1: Duality mapping in one dimension as a recombination 
of Majorana fermions. c(n) — [A{n) + iB{n)]/2 is the fermion 
defined on the original lattice site and d{n) = [B{n) +iA{n + 
f)]/2 is defined on the dual lattice site. 



it is inherently nonlocal. In the following, we generalize 
this duality transformation to two-dimensional systems 
and show that the transformation can be used to exactly 
map a classically ordered system to a topologically or- 
dered one. Two specific models are discussed. 

The first model is Wen's exactly soluble spin model 
defined on a square lattice^ 



H 
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where is the coordinate of lattice site. It is easy 

to see [Fij,Fi'ji] — and the model is thus exactly sol- 
uble. This model is shown to have robust topologically 
degenerate ground states and gapless edge excitations^. 

Let us first introduce the two-dimensional Jordan- 
Wigner transformation to fermionize the model^S 



j'<3 »' 



i' <i 



If we define the Majorana operators 



and By = i (4j - 



we find Fi, 



ij — ^ij^i+lj"-Bij_(-i_Bi_(-ij + i. 
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It is now interesting to generalize the duality to two 
dimensions and define fermions on vertical bonds (i, j) — 
(*,J + 1) 



4 = (Ay-zi?,,,+i)/2. 



It follows that 
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FIG. 2: (a) Graphical representation of Wen model. Each 
circle denotes one original lattice site on which two Majorana 
fermions A and B are defined. Each dotted rectangle con- 
tributes a ring-exchange term given by Eq. p2p . (b) After 
introducing fermions on vertical bonds, the model reduces to 
decoupled Ising chains. 



where fXij is related to the fermion d through a Jordan- 
Wigner transformation. The Wen Hamiltonian thus 
reads 



(16) 



This new Hamiltonian describes a set of decoupled Ising 
chains^. A first order phase transition from ferromag- 
netism to antiferromagnetism happens at g = 0. It is 
now straightforward to see that all two-point correlation 
functions are identically zero {cicr^) = in the ground 
state. This is so because CF\a\ contains Aij (or B^j+i) 
that is unpaired with its partner B^j+i (or A^) due to 
the fractionalization of aij into Aij and Bij and the re- 
combination of Aij and into /.tf^. 

It is also interesting to fermionize a Zeeman term 
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We notice that Ea. ([TC)) -l-Ea.pT |) is the same fermioinic 
Hamiltonian obtained by fermionizing quantum com- 
pass model using Jordan- Wigner transformation^*'. Af- 
ter including a Zeeman term, the Wen's soluble model 
is thus equivalent to the quantum compass model, which 
is shown to have dimensional reductio n^^'^'^ and a first- 
order phase transition at & = g^. 

The duality mapping can also be made explicit as fol- 
lows. Define fj-ij on the bond — + 1) 



\i'>i 




V 



(18) 



Let us first prove that /xf^ commutes with Without 
losing generality, let us assume I > j- Let us consider the 
overlaps between the original lattice sites involved in 
and /i^;. (a) If Z — j > 2, there is no overlap. Similarly, 
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FIG. 3: (a) Graphical representation of Ea. (|18p . /i^ is defined 
on a vertical bond as the product of spin components as la- 
belled on the origin lattice sites along the thick dark line, (b) 
Nearest neighbor coupling between fi^j and /if+i j. 

no overlap happens when I = j + I and k > i and two 
fi^ commute, (b) I = j + I and k ^ i, there is only 
one common site on which is involved in both fi^ . (c) 
I — j + 1 and k < i oi I = j (see Fig|3Kb), as far as 
a commutation relation is concerned, the only relevant 
part is 

(a^^al)^ [ayal), = Ka,^)^ • (19) 

Here, the subscripts 1 and 2 denote two different sites of 
original square lattice while r(red) and fc (black) are used 
to keep track of which /i^ the spin operators a come from. 
Based on the above consideration, we conclude that /i's 
are commuting with each other. It is also trivial to see 
(/Lt^) = /i^ and (/z^) = 1. Therefore, can be viewed 
as an Ising degree of freedom defined on a vertical bond. 
Let us now consider the interaction term between two 
nearest- neighboring as illustrated in FiglSjb), 

Therefore, the original Hamiltonian ^ is equivalent to 
Eq. (|16p and the Ising order is mapped to the quantum 
order studied by Wen^. The explicit mapping (fT8|) also 
allows us to examine the topological nature of the order- 
ing. In the basis of cr^, /if^ creates two kinks at the two 
ends and + 1) of the vertical bond. Therefore, 
the ordering is actually a condensation of kink-dipoles. 
The topological nature of the order in this model is thus 
similar to the one of Ising lattice gauge theory^^. 

The discussion of the exact mapping is so far limited 
to bulk terms. It is important to study the boundary 
conditions and demonstrate explicitly the dependence 
of ground degeneracy on the topology. In general, the 
boundary conditions will induce a nonlocal phase factor 
to the coupling strength between boundary spins, due 
to the phase term in the Jordan- Wigner transformation. 
The topology-dependent phase factor has profound con- 
sequence as we shall show shortly. One immediate conse- 
quence is that it determines the ground state degeneracy. 

The most interesting boundary condition is the case 
where we put the original spin model into a closed 



topology. A simple closed manifold is a torus by 
taking periodic boundary conditions along both direc- 
tions. The boundary term along y-direction is Hj^ = 
Ly'^i+i.Ly'^i+i I'^ti- It is clear that the phase term 
cancels and the periodic condition along the y-direction is 
mapped to a periodic boundary condition in the direction 
perpendicular to the Ising chain. The periodic bound- 
ary condition along x-direction induces a boundary term 
ffcrLj'^ij^ij-n'^Lj+i- This term is mapped 
to = gxJ2jt^tjt^L,,j '^it^ coupling strength 
given by 

5^ = 5 n '^Ij'^lj+i = 5 n ■ (21) 

The boundary term couples nearest neighboring chains 
nonlocally, which manifestly represents the hidden topo- 
logical structure in the original model. A direct conse- 
quence of this coupling is the partial lift of ground de- 
generacy. 

Another interesting case is the ribbon structure, where 
periodic boundary condition in y direction and open 
boundary condition in x direction are assumed. The 
open boundary condition in x direction is now mapped to 
the open boundary conditions in the spin chains. Conse- 
quently, the ground state now has an effect of dimension 
reduction and huge degeneracy 2^" , where we denote 
and Ly as the system sizes along x and y directions re- 
spectively. The degeneracy can actually be related to free 
Majorana fermions on boundaries. This can be shown 
by considering an equivalent geometry where we set pe- 
riodic boundary condition along x and open boundary 
condition along y. The ground state degeneracy is 2^^ 
in this case. The mapping thus leads to unpaired Ma- 
jorana fermions, Ai^L (j- = l,--,^a;) on the sites of top 
boundary and Bi^i(i = l,..,Lx) on the sites of bottom 
boundary. Ai^L^ are coupled to the bulk system through 

the boundary term in the form of Y[f=i ^i,Ly ■ Similarly, 
Bi i are coupled to bulk through Y[f=i-^i,i- Therefore, 
the operators that flip even umbers of unpaired Majo- 
rana fermions on top/bottom boundaries are conserved 
quantities and consequently lead to degenerated ground 
states. For instance, we can combine Ai^^ and ^2,l hito a 
fermion whose particle number (zAi^i^2,L + l)/2 is a con- 
served quantity. We thus have successfully mapped the 
global (nonlocal) Z2 degree of freedom of the decoupled 
Ising chains into a local degree of freedom of unpaired 
Majorana fermions at the ends of the chains. 

To illustrate our approach further, we show that sim- 
ilar physics can also be reached for the second Kitaev 
model defined on a honeycomb lattice-, 

H = - J2 E J^^J^k- (22) 

X—x.y,z X — bonds 

This topologically ordered model can be mapped to a 
model of spinless fermions whose ground states are char- 
acterized by local order parameters. 
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FIG. 4: Deformed honeycomb lattice and three types of 
bonds. 



Again, we fermionize this model using the Jordan- 
Wigner transformation. The idea is to deform the honey- 
comb lattice into a brickwall lattice as shown in FiglH We 
introduce the subscripts b and w to denote the white and 
black sites of a bond as illustrated in Fig[4l We also define 
the corresponding Majorana fermions = (c — c^)w/i 
and = {c + c^)^ for white sites, B}, = {c — c^)h/i and 
Ab = (c + c^)b for black sites. After a Jordan- Wigner 
transformation given by Eq.® and Eq. ((TO| . the Hamil- 
tonian of Kitaev model becomes 



H 



JxA^Afj ^ JyAf^A^ 

y — bonds 



X — bonds 



(23) 



z — bonds 



where a = iB\,Bw defined on each vertical bond. It is 
easy to see that a is a conserved quantity^^ and can now 
be taken as a number that can take either +1 or —1. 
The Hamiltonian is now quadratic in A and readily to be 
solved exactly. 

We are now ready to generalize the duality to brickwall 
lattice and introduce fermion on a z-bond, d = (A^, + 
i^h)/2 and S = (A„ — iAi,)/2, where A^j and A}, are 
the Majorana fermions on the white and black sites of a 
given z-bond. We thus have a model for fermions on a 
square lattice with site-dependent chemical potential 



4i? 



= JxJ2(d 



i 

-J.5]a.(24d,-1). 



i + Cy 



(24) 



Here By connects two ^-bonds and crosses a y-bond, sim- 
ilarly for Bx, as illustrated in FigHl This Hamiltonian 
describes a system of spinless fermions with p-wave BCS 
pairing and site-dependent chemical potential, where the 
ground states are characterized by local order parame- 
ters. Previous discussions about ground state degeneracy 
and vanishing two-point spin correlation functions can 
now be extended to this model straightforwardly. Un- 
paired free Majorana fermions also emerge naturally at 
open boundaries. For instance, a ribbon geometry can be 
achieved by breaking a row of z-bonds. For each broken 
z-bond, the Z2 degree of freedom a = iBj^Byj is fraction- 
alized into two unpaired free Majorana fermions, B}, — 
{c ~ c^)f,/i on the top boundary and i?^ = (c + c^)w on 
the bottom boundary. A detailed study of the fermion- 
ized Hamiltonian ([M)) will be presented elsewhere^^. 



In summary, we have successfully constructed exact 
mappings from topological orders to classical orders in 
two exactly soluble spin models. The topological de- 
pendence in the later model is manifestly represented in 
the terms resulted from the mapping of boundary condi- 
tions. Unpaired Majorana fermions on open boundaries 
and vanishing two-point spin correlation functions also 
follow naturally from our construction. Our work sug- 
gests a novel approach to construct certain topological 
orders from well-studied classical orders through a non- 
local transformation. Finally, we would like to point out 
that the transformation used in this work only works for 
a limited class of models. However, we conjecture that a 
general connection between topological orders and classi- 
cal orders might be possible and more beautiful nonlocal 
transformations are waiting to be discovered. 
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